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Several phenomena in high energy astrophysics have been recently related to clusters of galaxies and to cosmic 
ray interactions occurring inside these structures. In many of these phenomena the observable effects depend on 
the energy density of cosmic rays confined in the Intra Cluster (IC) medium, which is a poorly known quantity. 
We propose here that useful indications about this quantity can be obtained from present and future observations 
of galaxy clusters in the radio and hard X-ray frequency ranges. 



1. Introduction 

Clusters of galaxies are the largest gravita- 
tionally bound structures in the universe. They 
are exceptionally useful laboratories for cosmol- 
ogy and high energy astrophysics. 

From the cosmological side these structures 
probe: 

i) the amplitude and shape of the primordial fluc- 
tuation spectrum, because the mass (or tempera- 
ture) function of rich clusters is strongly sensitive 
to the value of the fluctuation power-spectrum, 
a{M,z). Thus it is possible to measure cr(M, z) 
on the cluster scales fitting the mass or temper- 
ature function to the local data (see, e.g., |l[] 
and references therein); 

ii) the evolution of baryons, both condensed in 
the form of galaxies and diffuse in the form of IC 
gas, which is abundantly present within the clus- 
ter potential wells. In fact, galaxies are respon- 
sible for the chemical enrichment of the IC gas 
at a level Fe/H ^ 0.3 of the solar value, during 
their life-cycles. Thus, high-quality X-ray spec- 
tral observations of galaxy clusters allow to study 
in details the physical state of the IC gas and the 
feedbacks from galaxy evolution ||] ; 

iii) the overall structure of the universe, using the 
possibility to measure directly Hq and flo through 
X-ray and Sunyaev-Zel'dovich effect (in the radio 



and sub-mm bands) cluster observations [Q. 

Clusters of galaxies are also relevant for high 
energy phenomena in large scale structures be- 
cause they can be regarded as the largest parti- 
cle accelerators in the sky. In the following, we 
will discuss some aspects in which galaxy clus- 
ters can yield important insights for high energy 
astrophysical phenomena. 

2. Clusters of galaxies as 7-ray sources 

There are different types and sites of particle 
acceleration in galaxy clusters from which 7-rays 
and neutrinos can be produced through the fol- 
lowing channels: 

p-\-p ^ tt" -\- X, 7r°^7 + 7 

P + P^TT^+X, TT^ ^ /i^l/p(P^), 

^i^ ^ i>f_,{Vf_,) + VeiVe) ■ (1) 

The possible sites are |P [H: i) normal galax- 
ies with low cosmic ray (CR) emitted powers, 
LcR ^ 10'*'^ erg/s; ii) active galaxies (like radio- 
galaxies living in the cluster cores or AGNs) with 
moderate powers, Lcr ~ 10''^ erg/s; iii) ac- 
cretion and/or merging shocks which are pro- 
duced during the cluster collapse and virializa- 
tion (see Fig.l) and can produce large powers, 
Lcr ~ 2-10''^ erg/s, through first order Fermi ac- 
celeration mechanism. In the pp collisions giving 



Figure 1. Shocks in large scale structures from N- 
body simulations [7] (courtesy of D.Ryu) in which 
are shown the distribution of the baryonic mat- 
ter (upper panel) and of the magnetic field (lower 
panel). Clusters of galaxies are found as the knots 
of the large scale matter condensations. 



rise to 7-ray production, bullet protons are sup- 
plied by the previous particle accelerators found 
within galaxy clusters, while the target protons 
are those of the hot (T ~ 10^ H- 10® K), diffuse 
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') and metal enriched IC gas (see 



g for a review) which is responsible for the clus- 
ter X-ray thermal emission. Thus, a correlation 
is expected between the X-ray (thermal) power: 
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and the 7-ray (non-thermal) power: 
Qi{E) oc YiappCnoQp{E)rl , 



(2) 
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where no is the central density of the IC gas, T 
its temperature and Tc the core radius of a cluster 
with density profile n{r) = no(l -I- x^)~^^^^, with 
X = r/rc and /3 ~ 0.6 0.8 (see [|| for details). 
The correlation among the X-ray and 7-ray fluxes 
(see Fig. 2) reads: 
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Figure 2. Gamma ray fluxes for clusters at z = 
0.023 in different cosmological scenarios: CDM 
(9.0 = 1; dashed curve), CDM+A (flo = 0.4; 
continuous curve), open CDM (ilo — 0.3; long- 
dashes curve) and mixed DM (flo = 1, i^iy = 0.3; 
dotted curve). The arrow indicates the ECRET 
upper limit for A 1656 ( Coma). 



In the previous formulae, Qp{E) is the CR 
spectrum which determines Lcr and conse- 
quently the cluster 7-ray emissivity, but it is a 
poorly known quantity in clusters of galaxies and 
in large scale structures in general. Assuming 
LcB. ~ lO'*^ erg/s (as produced by accretion 
shocks and/or active galaxies) we predicted that 
the contribution of galaxy clusters to the diffuse 
gamma ray background is 3%, quite indepen- 
dently of the assumed cosmological scenario and 
on cluster evolution (see Fig. 3). This result (ap- 
parently not relevant for the discussion on the 
origin of the DGRB) is indeed relevant because it 
shows that galaxy clusters cannot overcome the 
fraction of the DGRB left available from other 
(resolved or truly diffuse) sources summing up to 
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Figure 3. The expected DGRB from clusters, 
considering various sources of theoretical uncer- 
tainties in the cluster modelling. A flat CDM 
(T2o = l;h = 0.5; n = model is considered 
here. The thick solid line show the best fit to the 
EGRET data [9]. 



^ 80 -f- 95% (see discussion in g]). Larger frac- 
tions of the DGRB contributed by galaxy clus- 
ters would give sensitive constraints to the cur- 
rent scenarios for structure formation. 

Note, however, that these results rely on the 
(reasonable and motivated) value assumed for the 
CR power in clusters, Lcr ~ 10*'* erg/s, which 
depends on the CR spectrum Qp{E). So far, there 
are only heuristic derivations of such a crucial 
quantity. Nonetheless, in the following we will 
propose an operative procedure to obtain quanti- 
tative indications on the CR spectrum in galaxy 
clusters. 



3. Non-thermal emission from secondary 
electrons in galaxy clusters 

Most of the cosmic rays produced in galaxy 
clusters are confined in the IC medium for times 
larger than the age of the universe |^ How- 
ever, the maximum energy for which this con- 
finement is still efhcient depends strongly on the 
adopted diffusion model, so that for very extreme 
choices of the parameters, confinement can be 
limited to low energies. Here we shall consider 
the two limiting cases of complete confinement 
and no confinement in order to outline the gen- 
eral picture of the CR diffusion in clusters. Dur- 
ing the lifetime of the CR in the IC medium, they 
interact with the hot IC protons and produce sec- 
ondary particles through pp collisions (see eq.l). 
In such collisions charged pions are also produced 
and they produce a flux of neutrinos (see ^ [||) 
and electrons (here we will refer to both electrons 
and positrons as 'electrons', being this distinction 
irrelevant for our purposes). 

The global normalization of the gamma or neu- 
trino spectra obtained in previous works depends 
on a poorly known quantity, namely the energy 
density of CRs in clusters. Here we propose that 
useful indications on the CR energy density in 
clusters can be obtained by the radio measure- 
ments already existing and by X-ray measure- 
ments in the hard part of the spectrum, that 
will be available in the next future. In our 
model, radio and non-thermal X-ray emissions are 
produced by relativistic electrons through syn- 
chrotron and inverse compton scattering (ICS) 
off the microwave background photons, respec- 
tively. 

We assume here a CR spectrum of the form: 
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where Eq ^ 1 GeV and Lp is the cluster total CR 
luminosity. The rate of production of secondary 
electrons from pp collisions is pcJ : 
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in the case of complete CR confinement. For the 
case of no confinement we obtain: 



,XE,) = 36 



LpUHcao — - 



^(7g,^?)(7g-l)(7g-2) 



6Do ml 



-E. 



-Ig-V 



(7) 



Here nn is the average gas density in clusters, 
(To = 3.2 X fO-26cTO2, to - '^"'^ ^ ini7^-i 



2.06 X W^h~- s 



is the age of universe and ra^^ and are 
the pion and muon masses. The diffusion co- 
efficient can been written as D{E) = DqE^. 
The functions /(7g) and I{^g^rj) are defined 
as /(7g) = /g^ (ix/^(a;)x'^9"2 and I{jg,v) = 

dxUix)x'^o+'^-^, where U{x) = f.34(f - 
x)^'^ + e~^^^, is a scaling function entering the 
calculation of the secondary electrons from the 
decay of charged pions (see [0). In the follow- 
ing we use 7g = 2.4, checking at the end how the 
results change for jg — 2.1. 

The diffusion coefficient enters effectively in the 
calculation of the case in which no confinement is 
realized. For the case of a Kolmogorov spectrum 
of fluctuations, it can be easily written in terms 
of the average magnetic field B, the coherence 
scale Ic and the CR energy E (see 0): 



D{E) 



Ir. 



3 X io-i.V35;V3 ^^^^^^^^ 



2/3 



CTO^ (8) 



where E is the CR energy in GeV and is the 
magnetic field in fiG. The coherence scale, Ic, of 
the magnetic field is a quite uncertain parame- 
ter: we assumed for it an upper value, 430 kpc, 
corresponding to the case in which the magnetic 
field lines are stretched by the bulk motion of the 
galaxies in the clusters on a typical scale equal to 
the average galaxy separation. 

At the electron energies we are interested in, 
the time for energy losses is smaller than the typ- 
ical diffusion time, so that the equilibrium elec- 
tron spectrum is ne(-E'e) Qe{Ee)Te{Ee), where 
Te{Ee) = Ee/{dEe/dt) is the time scale for losses. 
The details of the calculation of the radio and X- 
ray emission from synchrotron and ICS are given 



in ; here we shall limit ourselves to outline the 
results obtained in different cases. 

In the case of complete CR confinement, the 
radio power spectrum reads: 



1.7 X lO^'^Li^n^ih-^viHzy 



erg s 



-1 Hz-\ (9) 



with B = 1 fiG and 7g = 2.4. Here L44 is the 
cluster CR luminosity in units of lO**^ erg/s, 713 
is the gas density in units of 10^"^ cm^'^ and ft.75 is 
the Hubble constant in units of 75 km/{s Mpc). 
With the same parameter choice, the X-ray spec- 
trum reads: 

IxiEx) = 

5.3 X 10^°LMn3hj^^ExikeV)'^-'^s~^ fce^-^ (10) 

While the X-ray flux is very weakly dependent on 
the magnetic field, the radio flux changes strongly 
with B. If we take B « O.I ^iG eq.(|) yields; 

1.2 X lO^^LMnahj^^viHz)-^-^ erg 5"^ fl"z-^(Il) 

This fact has a relevant implication: an accurate 
measurement of the flux in the hard X-ray region 
of the cluster spectra can fix the product -^44713 
so that, in turn, a measure of P, (j^) can give in- 
formations on the value of B. 

In the case of no confinement we obtain for the 
radio power spectrum: 

5.6xlO^'^Uin3h^liy{Hz)-^-^^^ erg s'^ Hz-\{12) 
and for the X-ray spectrum: 

IxiEx) = 

l.AxlO^°Liin3h^lEx{keV)-^-^^^s^^ A:ey-^(13) 

Again Ix{Ex) does not change appreciably with 
the magnetic field, while Pr{iy) does. For B = 
0.1 iiG we obtain: 

l.SxlO^^LiiUsh^^i^iHz)-^-^^^ ergs^^ Hz-\{U) 
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4. A future outlook 

At present, there are only a few observations 
of cluster radio halo spectra, the most detailed 
one referring to the Coma cluster Based 
on the available observations, we can only put 
weak constraints to the cluster CR energy den- 
sity, uJcB. ~ 10^^ eV/cm^ (confinement with B = 
l^G) or uj^ji ^ 1 eV/cm^ (confinement with B = 
O.l/xG). Moreover, note that the EGRET upper 
limit on the 7-ray emission from Coma (see Fig. 2) 
yields ^0.1 eV/cm^. 

From this first comparison of our results with 
the available data, it is evident the extreme im- 
portance of a definite detection of non thermal 
X-ray tail in the spectra of galaxy clusters. In 
fact, such X-ray spectra are very weakly depen- 
dent on the value of the magnetic field, so that 
X-rays can in principle fix the value of L^^ns, pro- 
vided that the radio halos observed in the same 
clusters are due to the same secondary electrons. 
Then, in turn, radio observations can be used to 
give reliable indications concerning both the am- 
plitude and the power spectrum of the magnetic 
field in clusters. 

The available observations at £' > 20 keV 
from Coma can only provide upper limits to the 
hard, non-thermal X-ray tails . In our model, 
these limits yield ^44713 ^ 60 ^ 600 where the 
two values refer to the extreme cases of confine- 
ment and no confinement, respectively. These 
give L44 1.4 • 10^ ~ lA ■ 10^ for an average IC 
gas density in Coma of 77.3 ~ 0.43. The previ- 
ous limits are very weak and do not allow, at the 
moment, neither to distinguish between the two 
regimes nor to strongly constrain the CR energy 
density in Coma. The results obtained for the 
case of an electron production spectrum with in- 
dex 7g = 2.4, do not change appreciably if we use 
7g = 2.1, which is somewhat close to a lower limit 
for this slope if the parent CRs are assumed to be 
produced by acceleration at supernovae shocks. 

Beyond the details of the calculations we pre- 
sented here, the possible detection of hard X-ray 
tails in the spectra of galaxy clusters will provide 
new insights for the role of high-E phenomena 
in large scale structures. This is at hand in the 
hard X-ray region with the data on galaxy clus- 



ters achievable with the SAX satellite, which is 
now operating. In fact, the PDS instrument on 
board SAX has a limiting flux ^ 1Q~^^ erg /crri^s, 
with AE ~ 15% (at 60 keV) and a broad band, 
AE ^ 15 300 keV, response which makes it 
possible to detect hard X-ray tails in the spectra 
of several nearby clusters like Coma, A2163 and 
A2199. 

Besides the analysis of the two extreme cases of 
CR confinement we presented here, the detection 
of a possible steepening in the X-ray tail spectra 
can suggest (see sect. 3) that CR begin to be un- 
confined above some threshold energy, while no 
effect should be present in the radio spectrum. 
A more detailed discussion of this effect will be 
reported in a forthcoming paper | pl] |. 

Thus, in the light of the present results, better 
constraints on the CR energy density in the IC 
medium are expected to become available in the 
next future, together with a better insight on the 
structure of the IC magnetic fields. 
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